Trapping of Neutral Mercury Atoms and Prospects for Optical Lattice Clocks 
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We report a vapor-cell magneto-optical trapping of Hg isotopes on the ^5*0 — ^Pi intercombination 
transition. Six abundant isotopes, including four bosons and two fermions, were trapped. Hg is the 
heaviest non-radioactive atom trapped so far, which enables sensitive atomic searches for "new 
physics" beyond the standard model. We propose an accurate optical lattice clock based on Hg 
and evaluate its systematic accuracy to be better than 10~^^. Highly accurate and stable Hg-based 
clocks will provide a new avenue for the research of optical lattice clocks and the time variation of 
the fine- structure constant. 
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Precision tests of fundamental symmetries with atoms 
and molecules play an important role in probing new 
physics beyond the standard model of elementary par- 
ticles. The "new physics" sensitivity of such experi- 
ments depends steeply on the nuclear charge Z of the 
atom. In particular, parity- violating amplitudes change 
as p,], atomic CP- violating (T-, P-odd) permanent 
electric dipole moments (EDM), as Z^ [2], and CP- 
violating polarizabilities, as Z^ [3]. As a result, the most 
accurate limit to date on the atomic EDM has been set 
with Hg atoms {Z = 80) and on the electron EDM with 
Tl {Z = 81) atoms [4, 5]. In a quest for more sensitive 
probes, experiments with heavier atoms, Fr {Z = 87) and 
Ra {Z = 88), have been proposed [6, 7]. However, all the 
elements with Z > 83 are radioactive and consequently 
experience associated drawbacks. 

In this Letter, we report the vapor-cell magneto-optical 
trapping (MOT) of Hg isotopes and discuss its applica- 
tion to accurate optical lattice clocks [8]. Located at 
the very edge of the stability valley, Hg is, so far, the 
heaviest non-radioactive atom that has been laser-cooled 
and trapped. It is anticipated that these ultracold heavy 
atoms will enable experimental searches of high-Z atom 
physics with unprecedented sensitivity due to reduced 
motional effects, longer interrogation time, and higher 
optical thickness. Moreover, as demonstrated for Yb iso- 
topes [sl , a remarkable variety of trapped Hg isotopes (4 
bosonic and 2 fermionic ones) may allow attaining Bose 
or Fermi-degeneracy, which would further improve the 
measurements' sensitivity. 

We systematically evaluate various sources of uncer- 
tainty for the Hg-based optical lattice clock and argue 



that an accuracy of better than 10~^^ is attainable, which 
is an order of magnitude of improvement over Sr [sl or 
Yb-based clocks because of reduced susceptibility to 
the blackbody radiation field, which sets a niajor limi- 
tation on the accuracy of atomic clocks [13, EH. This 
projected accuracy is competitive with that of the best 
ion clock with A1+ [11]. In particular, we take advan- 
tage of the high nuclear charge Z of Hg to explore the 
variation of the fine-structure constant a that generically 
controls the strength of relativistic effects in the atom as 
(aZ)^. While astrophysical observations suggest that a 
may drift at the rate of a/a ~ 10~^^yr~^ [12], recent 
laboratory tests based on the frequency ratio of accu- 
rate ion clocks and Cs clocks set the upper bound of 
\a/a\ < 1.3 X IQ-^^yr-i limited by Cs clock' s mac- 
curacy of 1 X 10~^^. Significant enhancement in measure- 
ment precision will be expected for direct comparison of 
optical clocks with frequency comb technique Q: the 
frequency change Su of Hg-based clocks operated at uq 
can be measured in reference to Sr [15] or Hg+ [13] clocks, 
where Hg, Sr, and Hg+ have sensitivity {5u/uo)/{5a/a) 
of 0.81, 0.06, and —3.19, respectively, to the fractional 
change of Sa/a 0, 0|. 

The laser cooling and clockwork with neutral Hg ex- 
ploit its divalent electronic structure and parallel tech- 
niques developed for Sr [1, [l^. The energy levels of 
Hg are shown in Fig. 1. The spin- forbidden ^5*0 — ^Pi 
transition at 253.7 nm with the natural linewidth of 
7/27r = 1.3 MHz [Til provides an adequate cooling tran- 
sition down to the Doppler temperature of 31 /iK. The 
^So — ^Pq transition at 265.6 nm, which is weakly allowed 
for odd isotopes with a nuclear spin (/ 7^ 0) or even iso- 
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FIG. 1: Energy- level diagram for Hg. The ^5*0 — ^Pi transi- 
tion at 253.7 nm is used for magneto-optical trapping of the 
atoms. The ^5*0 — ^Po transition at 265.6 nm serves as a clock 
oscillator. The "magic wavelength," Al = 360 nm, is the far- 
off resonance for the relevant transitions originating from the 
^So and ^Po levels. 



topes (/ = 0) dressed by external fields fio', "2^, will be 
exploited as a clock transition. The "magic wavelength," 
where the light-shift perturbations of the clock levels are 
canceled out [1, [2l|, is calculated to be Al = 360 nm, as 
discussed later. 

Trapping neutral Hg is experimentally challenging due 
to limited laser power available at deep-UV cooling tran- 
sition wavelength. To overcome this difficulty and take 
advantage of its heavy mass and the large radiation pres- 
sure at the shorter wavelengths, we used a vapor-cell 
MOT [22] to capture the atoms from the low- velocity tail 
of the Maxwell-Boltzmann distribution. Our apparatus 
is shown in Fig. 2. A 253.7 nm light was generated by 
frequency quadrupling a laser diode system operated at 
1014.9 nm. A master laser was an external-cavity laser 
diode locked to a reference cavity to reduce its linewidth 
to less than 100 kHz. This master laser was amplified by 
a tapered amplifier up to 1 W and injected into a bow- 
tie cavity with a PPLN crystal to generate 240 mW at 
507.4 nm. This was further frequency-doubled by a BBO 
crystal to obtain roughly 10 mW at 253.7 nm. 

This UV light was fed into a vacuum chamber with its 
base pressure of 3 x 10~^^ Torr to form a vapor-cell MOT 
with a folded beam configuration with a single stroke. 
The 1/e^ diameter of the trapping beam was 20 mm with 
a power density of 0.3/o, where Iq = lOmW/cm^ is the 
saturation intensity of the ^5*0 — ^Pi transition. After 
introducing Hg, the vacuum pressure was increased to 
8 X 10~^ Torr that indicated the Hg-atom number density 
of 7 X 10^ cm~^. The typical magnetic field gradient of 
a quadrupole magnetic field was 10 G/cm along its axis 
or the Zeeman shift gradient of 21 MHz /cm for the ^Pi 
state. 

Hg has 6 stable isotopes with abundances of more than 
6%; 4 bosonic (^^^Hg, ^ooHg, '^^'^Rg, and ^o^Hg) isotopes 
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FIG. 2: Experimental setup. A 253.7 nm UV light was gen- 
erated by frequency quadrupling a diode laser operated at 
1014.9 nm with nearly 1 W. A 10 mW UV light was fed into 
a vacuum chamber to form a vapor-cell MOT with a folded 
beam configuration. ECLD: External Cavity Laser Diode, 01: 
Optical Isolator, TA: Tapered Amplifier, PZT: Piezo-electric 
Transducer, QWP: Quarter Wave Plate. 



with nuclear spin 7 = and 2 fermionic ^^^Hg and ^^^Hg 
isotopes with 7 = 1/2 and 3/2, respectively. The spectra 
of trapped Hg isotopes are shown in Fig. 3, where the 
UV frequency was swept toward higher frequencies with 
0.5 MHz/s. The typical atom filling time was measured 
to be ~ 1 s for ^^^Hg. The trap spectra are about 10 MHz 
or 87 wide, whose asymmetry are explained in terms of 
atom loading dynamics into the MOT [22, 23]. While 
atom capturing efficiency is higher in the red wing, the 
trapped atoms scatter more photons as the laser detuning 
A becomes smaller. The atoms suddenly boil up from the 
MOT when they cross the Doppler cooling condition of 
A = —7/2 over to their resonance, corresponding to the 
sharp cutoff in the blue wing. 

In Fig. 3, the number of trapped isotopes was nearly 
proportional to their natural abundance, indicating that 
atoms can be further accumulated until two-body light- 
assisted collisions limit their density. For odd isotopes 
(with a hyperfine structure), stable MOTs were observed 
for the ^So{F = 1/2) ^ ^Pi{F = 3/2) transition for 
i^^Hg and ^So{F = 3/2) ^ ^Pi{F = 5/2) for ^oiRg, 
as expected [23]. From the fluorescence of the MOT, the 
number of trapped atoms was roughly 10^ for ^^^Hg with 
a 1/e^ cloud radius of about 0.2 mm. The atomic tem- 
perature was measured to be less than 50 jaK by the time- 
of- flight technique, which was consistent with its Doppler 
temperature. Even lower temperatures are expected for 
odd isotopes because of the Sisyphus cooling mechanism 
^24] . Since atom loading efficiency into the vapor-cell 
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MOT strongly depends on trapping laser diameter [22[, 
substantial increase in the number of trapped atoms may 
be expected by applying brighter light sources [25[. 
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FIG. 3: Fluorescence spectrum of magneto-opt ically trapped 
Hg isotopes as a function of the laser frequency at 253.7 nm. 
The bars at the bottom indicate the positions of the reso- 
nances of individual isotopes. Detailed profiles are drawn us- 
ing the magnified frequency scale shown at the top. 



In the following, we discuss an ultrastable optical lat- 
tice clock operating on the 6s^ ^Sq — 6s6p ^Pq transition 
(see Fig. 1) of Hg. Starting with the magneto-optically 
trapped atomic samples, as described above, a possi- 
ble experimental scenario is as follows: a 100-/iK-deep 
ID (one-dimensional) optical lattice is prepared inside 
a power build-up cavity with a peak intensity of the 
standing wave of = 1 MW/cm^, where the light- 
shift Sul is given by Sul/Il ^ 2 kHz/(kW/cm^). The 
trapped atoms are then sideband-cooled to the vibra- 
tional ground state of the fast axis on the ^5*0 — ^Po 
clock transition with an appropriate relaxation process 
via the state (see Fig. 1). The optical lattice clock 
can be operated after adiabatically reducing the lattice 
intensity down to II = 50 kW/cm^, which corresponds 
to the potential depth of ~ 10 En, in order to minimize 
the hyperpolarizability effects of the lattice lasers. Here, 
Er = {h/ Xl)'^ /2m is the photon recoil energy of Hg 
atoms scaled by the magic wavelength A^. In Ref. [26], 
it is pointed out that optical lattices as shallow as 10 En 
could well meet the 10"-^^ accuracy goal by suppressing 
the tunneling of atoms with the gravitational acceleration 
when the ID lattice is vertically oriented. 

We evaluate the accuracy of the Hg optical lattice clock 
by carrying out relativistic calculations of the relevant 
atomic properties. Hg has 80 electrons, and the correla- 
tions play an essential role in characterizing the atomic 
properties. We employ the method of configuration- 
interaction (CI) coupled with the many-body perturba- 
tion theory (MBPT) \2?\. The basis set was constructed 
in the V*^^"^) approximation. In CI+MBPT, the posi- 



tions of the low-lying energy levels are reproduced within 
a few 0.1%. Hyperfine-structure constants and dipole 
matrix elements have a theoretical accuracy better than 
10%. The most relevant test for the prediction of the 
"magic wavelength" is the difference in the static polar- 
izabilities of the 6s6p ^Pi and the ground states, measured 
in Ref. [28] to be 26.66(47) a.u. Our computed value of 
26.95 a.u. is in agreement with this experiment. 

We start by determining the natural width of the clock 
transition. Although a single-photon 6s^ ^Sq — 6s6p ^Pq 
transition is forbidden by the conservation of the an- 
gular momentum, hyperfine coupling with the nuclear 
angular momentum / 7^ opens a weak decay chan- 
nel [29]. Our calculations of the hyperfine quenching re- 
sult in the rate coefficients A{^^^I{g) = 1.3 x 10~^s~^ 
and A(^^^Rg) = 8.8 x 10"^ s-^ In even isotopes (/=0), 
the natural linewidth would be introduced only due to 
highly suppressed El-Ml two-photon decays. 

In an optical lattice, the shift of atomic levels is propor- 
tional to the a.c. polarizability a (A), and we compute it 
for both clock levels. They cross at the magic wavelength 
of Xl ^ 360 nm, which agrees with the semi-empirical es- 
timates of 358 nm [l5[ and of 342 nm [30] . Compared to 
the El contribution, the higher-rank multipolar a.c. po- 
larizabilities are suppressed, and they merely shift by 
a small fraction. 

In addition to the dominant scalar contribution, the 
fermionic isotopes (/ 7^ 0) acquire vector ^'^"^(Al) (and 
tensor for /> 1/2) components of the a.c. polarizability. 
The non-scalar contributions are suppressed by the ratio 
of the hyperfine interaction to atomic energy. The ^Pq 
level is affected the most due to the nearby levels of the 
^Pj fine-structure manifold. We find ^'^''^(Al) = 0.040 
a.u. for i^^Hg and ^^^^(Al) = -0.045 a.u. for ^oiHg. 
Such vector shifts can be eliminated using an averaging 
procedure [3l| : Since the vector light shift is proportional 
to the of the Zeeman sublevels, alternating measure- 
ments of f± of the ^So{F, :^mF)—^Po{F, ±mi?) transition 
frequencies and taking the average /o = (/+ + /_)/2 can 
significantly reduce the fractional contribution p of the 
vector light-shift down to Suy ^ I'xp^ 1 where Vx ^ 50 kHz 
is the lattice trapping frequency for the 10£^i^-deep lat- 
tice. By applying nearly linearly polarized light with its 
ellipticity of less than 10~^, 5vy ~ 10~^ Hz makes a neg- 
ligible contribution to our goal accuracy. 

We also made a semi-empirical estimate of the second- 
order (in the lattice intensity II) corrections to the 
a.c. Stark shift of the clock frequency 5vh = 3?(/3)/2, 
where (5 = (71^0 — 73Pq)/64 is the difference between 
the clock-level hyperpolarizabilities and 3?(/3) denotes its 
real part. The calculations in the single-electron ap- 
proximation, as carried out in Sr [8|, infer that 73 
is about 10^ times larger than 71 5^ at A^ = 360 nm. 
It becomes smaller with decreasing A^ and crosses zero 
near 347 nm for a linear polarization. Numerically, 
(3 = -(103 + 5.5z) mHz/(MW/cm^)2 for Al = 360 nm. 
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and f3 = -(21.7 + 1.5i) mHz/(MW/cm^)2 for Xl = 351 MinES. 
nm. The imaginary part of P determines the two-photon 
ionization rate P2ph = — 47r^(/3) /£. For Il= 50 kW/cm^ 
and the linear polarization, the induced shift 5uh is below 
0.3 mHz, and the two-photon ionization rate does not ex- 



ceed 2 X 10 ~^ s~^, which is smaller than the off-resonant 
(single-) photon scattering rate of 5 x 10~^ s~^ of the lat- 
tice lasers. For the circular polarization, both the real 
and the imaginary parts of P are several times larger, 
without any node in the indicated region of wavelengths. 
These second-order corrections are sensitive to both cor- 
relations and the magic wavelength of A^; therefore, they 
ultimately need to be investigated experimentally. 

The sensitivity to magnetic fields B is determined by 
the difference in the Lande g-factors of the ground and 
excited states. We find ^^(-^^^Hg) = —4.7 x 10~^ and 
dgC^^^Rg) = 1.8 X 10-"^ with the CI+MBPT calculations. 
While the linear Zeeman shift can be eliminated by the 
same averaging protocol as described previously [3l| , the 
simultaneously obtained Zeeman shift A/ = /+ — /_ with 
Sg given above provides a real-time probe of the applied 
magnetic field to correct the second-order Zeeman effects 
of Siy^z^ = -(3zB^ with (3z = 24.4mHz/Gl 

Finally, at room temperature, the accuracy of the 
clock is affected by blackbody radiation (BBR). The in- 
duced shift of the clock transition frequency, Scobbr ^ 
— j^{a7T)^T^ X (ofspQ — Q^i^o)? depends on the difference 
in the static polarizabilities of the clock states. Our 
CI+MBPT calculations result in asp^ = 54.6 a.u. and 
aisQ = 33.6 a.u., leading to the BBR shift of ^^^bbr ^ 
-0.181 Hz at T = 300 K. It is noteworthy that the frac- 
tional correction (^^^bbr/^o ^ — 1.60 x 10~^^ is an order of 
magnitude smaller than the Sr and Yb based clocks [l^ . 
Assuming a temperature uncertainty of 0.1 K, the frac- 
tional uncertainty due to the BBR is 2 x 10~^^. 

In conclusion, we have demonstrated a magneto- 
optical trapping of six isotopes of the Hg atom, the heav- 
iest non-radioactive atom trapped so far. We have shown 
that Hg is a promising candidate for highly accurate op- 
tical lattice clocks with an estimated inaccuracy of less 
than 10~^^, which is mainly limited by the hyperpolar- 
izability effects, higher-rank multipolar polarizabilities, 
and BBR shifts. Although, current performances of op- 
tical lattice clocks have been investigated at the 10 

level [m, m, [13 limited by the uncertainty of Cs refer- 
ence clocks, a highly accurate and stable Hg-based clock 
will provide a new avenue for improved systematic studies 
of optical lattice clocks, such as the precise determination 
of the "magic wavelength" and other properties beyond 
the 1 X 10~^^ uncertainty level. Furthermore, because of 
its high Z, Hg may allow a sensitive search of the tem- 
poral variation of a combined with other optical clocks, 
e.g., Sr-based clocks as a possible anchor [15']. 
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